Torsion and accelerating expansion of the universe in quadratic gravitation 
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Several exact cosmological solutions of a metric-affine theory of gravity with two torsion functions 
are presented. These solutions give a essentially different explanation from the one in most of 
previous works to the cause of the accelerating cosmological expansion and the origin of the torsion 
of the spacetime. These solutions can be divided into two classes. The solutions in the first class 
O ; define the critical points of a dynamical system representing an asymptotically stable de Sitter 

spacetime. The solutions in the second class have exact analytic expressions which have never 
been found in the literature. The acceleration equation of the universe in general relativity is 



only a special case of them. These solutions indicate that even in vacuum the spacetime can be 
endowed with torsion, which means that the torsion of the spacetime has an intrinsic nature and a 



■ geometric origin. In these solutions the acceleration of the cosmological expansion is due to either 
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the scalar torsion or the pseudoscalar torsion function. Neither a cosmological constant nor dark 

5— i ' 

energy is needed. It is the torsion of the spacetime that causes the accelerating expansion of the 
universe in vacuum. All the effects of the inflation, the acceleration and the phase transformation 
from deceleration to acceleration can be explained by these solutions. Furthermore, the energy 

> . 

, and pressure of the matter without spin can produce the torsion of the spacetime and make the 

' . ^ , expansion of the universe decelerate as well as accelerate. 
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, I. Introduction 

In the last few years the realization that the universe is currently undergoing an accelerated expansion phase and 
the quest for the nature of dark energy has renewed interest in so-called modified gravity theories (for a review see 
[1]). In these theories one modifies the laws of gravity so that a late-time accelerated expansion is produced without 
recourse to a dark energy component, a fact which renders these models very attractive. The simplest family of 
modified gravity theories is obtained by replacing the Ricci scalar R in the usual Hilbert- Einstein Lagrangian with 
some function f(R) (for reviews, see , [2-5]). 

There are actually three versions of f(R) gravity: Metric f(R) gravity, Palatini f(R) gravity, and metric-affine 
f(R) gravity. In fact, these are physically different theories rather than manifestations of the same theory in different 
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guises, as the different variational principles yield inequivalent equations of motion (except when the action is the 
Einstein-Hilbert and matter is minimally coupled to geometry). In metric f(R) gravity, the action is varied with 
respect to the metric as usual (for an introduction see [6]). Palatini f(R) gravity comes about from the same action 
if we decide to treat the connection as an independent quantity. The connection, however, does not enter the matter 
action. Such a approach were introduced and initially studied by Buchdahl [7] and has attracted a lot of interest as 
possible infrared modifications of general relativity (for a shorter review of metric and Palatini f(R) gravity see [8]). 
It has recently been generalized to f(R) theories with non-symmetric connections, i.e. theories that allow for torsion 
[9] and f(R : R tlu R ilv ) theories [10]. In mctric-affine f(R) gravity the matter action is allowed to depend also on the 
connection. In addition, the connection can include both torsion and non-metricity [11]. 

It has been shown that even in the most general case of Palatini /(i?) gravity where both torsion and non-metricity 
are allowed, the connection can still be algebraically eliminated in favor of the metric and the matter fields [12]. 
Clearly, f(R) actions do not carry enough dynamics to support an independent connection which carries dynamical 
degrees of freedom. However, this is not a generic property of generalized Palatini gravity. The addition of the 
RfivR^ term to the Lagrangian radically changes the situation and excites new degrees of freedom in the connection. 
The connection (or parts of it) becomes dynamical and so, it cannot be eliminated algebraically. If the connection is 
torsion free, the dynamical degrees of freedom reside in the symmetric part of the connection [13]. 

In generic metric-affine theories the addition of the R^R^ term to the Lagrangian makes the propagating degrees 
of freedom reside in both the antisymmetric and symmetric parts of the connection. In other words, the dynamical 
degrees of freedom can be both torsion and non-metricity. In these theories torsion field plays a fundamental role: it 
contributes, together with curvature degrees of freedom, to the dynamics. Propagating torsion is the key feature of 
these theories [14, 15]. 

Torsion proves to be essential for total angular momentum conservation when intrinsic spin angular momentum is 
relevant (for reviews on torsion, see [16, 17, 18]). It has been argued that torsion must be present in a fundamental 
theory of gravity [19, 20]. In the teleparallel gravity, for example, torsion plays a central role (for a shorter review 
see [21]). Recently, models based on modified teleparallel gravity, namely f(T), were presented. In these models the 
torsion proves to be the responsible of the observed acceleration of the universe [22]. 

The rediscovery of the metric-affine (Palatini) formulation was mainly driven by the interest in finding cosmological 
scenarios able to explain the current observations. Using the respective dynamically equivalent scalar-tensor represen- 
tation of Palatini f(R) gravity some cosmological models with asymptotically de Sitter behavior have been presented 
[23] . It was shown that adopting the metric-affine formulation together with an action that includes a term inversely 
proportional to the scalar curvature, such as the one in [24], can address the problem of the current accelerated 
expansion equally well as when using the purely metric formalism [25]. Additionally, it was found that f(R) theories 
of gravity in the metric-affine formulation do not suffer from the problems for the metric formulation. On the other 
hand, although cosmology in the theories with the Lagrangian including R 2 and R^ V R^ terms have been studied in 
the purely metric formulation (for example see [26]) and the Palatini formulation [27], the similar cosmological models 
in the metric-affine formulation have not been discussed thoroughly in the literature. Especially, the cosmological 
effect of torsion in metric-affine theories of gravity has not been explored extensively. We have not known whether 
the dynamical torsion could lead to a de Sitter solution and then be used to explain the observed acceleration of the 
universe. An answer will be given in this paper. 
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The metric-affine approach has been widely used in order to interpret gravity as a gauge theory many times over 
the years (see, for example, [28] for a study on f(R) actions and [29] for a thorough review). In recent years it has 
bee used in cosmology to interpret the accelerating expansion of the universe [30, 31]. In this approach the structure 
of the gravitational equations and physical consequences of cosmology, in particular, the situation concerning the 
accelerating expansion depend essentially on the form of the Lagrangian. The metric-affine gravity can be divided 
into different sectors in dependence on the number of nonvanishing components of the torsion tensor and the order of 
the differential equations. One sector of the metric-affine gravity is so-called dynamical scalar torsion sector considered 
in [30] . Starting from a Lagrangian consisting of R 2 and the quadratic torsion terms a cosmological model has been 
constructed. This model can contribute an oscillating aspect to the expansion rate of the universe. A different model 
of acceleration with torsion but without dark matter and dark energy has been presented in [31]. The Lagrangian of 
it is the most general form including the linear in the scalar curvature term as well as 9 quadratic terms (6 invariants 
of the curvature tensor and 3 invariants of the torsion tensor with indefinite parameters). Its Lagrangian involves 
too many terms and indefinite parameters, which make the field equations complicated and difficult to solve and the 
role of each term obscure. In order to simplify the field equations some restrictions on indefinite parameters have to 
be imposed. Under these restrictions, especially, all the higher derivatives of the scale factor are excluded from the 
cosmological equations. The question is whether such a complicated Lagrangian is necessary. Can we use a simpler 
Lagrangian to construct a model of cosmic acceleration? In fact all the indefinite parameters in the Lagrangian in [31] 
have been combined into four new ones, which implies that some terms are not necessary and the Lagrangian can be 
simplified. In this paper we will show that a rather simpler Lagrangian, R + aR 2 + PR flv R llv +"fT tl l/p T ' ^ p , is sufficient 
and necessary to construct a model of cosmic acceleration. The terms pR^ v R^ v and jT p vp T p p play different roles 
in the theory: the former determines the structure of the field equations while the latter determines the behavior 
and the stability of the solutions. The f3R pu R pv term leads to different structure of the cosmological equations from 
the one in [30]. In addition to the simplicity the main advantage of this Lagrangian is to permit exact or analytic 
solutions which have not been found in previous works. For any physical theories, to find exact or analytic solutions is 
an important topic. Next comes the physical interpretation of the solutions thus obtained. Mathematically de Sitter 
spacetimc as the maximally space is undoubtedly important for any gravity theories. From the observational side, 
recent studies illuminate that both the early universe (inflation) and the late-time universe (cosmic acceleration) can 
be regarded as fluctuations on a de Sitter background. So de Sitter solutions take a pivotal status in gravitational 
theories, especially in modern cosmology. 

We will follow the approach of [26, 30,31] rather than the one in [27] to avoid getting involved in debate on the 
transformation from one frame to anothcr[32]. We choose the tetrade^ and the spin connection T IJ p instead of the 
metric g pu and the affine connection r A pv as the dynamical variables following the gauge theory approach [30]. The 
descriptions in terms of the variables (e 1 pi T JK U ) and (<7 M „, T x pa ) are equivalent in our approach (the argument in 
detail see [16]). We will concentrate on the role of torsion subject to the metricity. In this case only the torsion part 
of the connection is independent of the metric (or tetrad). 

Because the field equations can result of order higher than second and very difficult to handle, the theory of 
dynamical systems provides a powerful scheme for investigating the physical behavior of such theories [33] for a wide 
class of cosmological models. The dynamical system approach has acquired great importance in the investigation 
on various theories of gravity. Some works have been done in the case of scalar fields in cosmology and for scalar- 
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tensor theories of gravity [34]. This approach has the advantage of offering a relatively simple method to obtain 
exact solutions (even if these only represent the asymptotic behavior) and to obtain a (qualitative) description of 
the global dynamics of the models. Such results are very difficult to obtain by other methods. The application of 
this method has allowed new insights on higher order cosmological models and has shown a deep connection between 
these theories and the cosmic acceleration phenomenon. It makes possible not only to develop experimental test for 
alternative gravity but also to allow a better understanding of the reasons underlying the success of the theory. The 
dynamical systems approach has been used to investigate universes in theories of gravity [26,30, 31]. In contrast 
with [31] we allow the field equations to contain higher derivatives. We will see that for the Lagrangian of the form 
R + aR 2 + PR I1V R IXV + 7T M vp T p vp the field equations can be simplified and solved exactly for some choices of a, f3 
and 7. Some meaningful consequences can be inferred from the solutions obtained. The accelerating expansion of the 
universe can be explained without a cosmological constant or dark energy. A vacuum spacetime can possess torsion 
which causes the acceleration of the cosmological expansion. The conception of vacuum as physical notion is changed 
essentially. Instead of it as passive receptacle of physical objects and processes, the vacuum assumes a dynamical 
properties as a gravitating object. The torsion of the spacetime can be produced by the energy and pressure besides 
the spin of matter. 

The paper is organized as follows. In section II the gravitational field equations are derived following the approach of 
[29, 30, 31]. Using them to the spatially flat Friedmann- Robertson- Walker metric a system of cosmological equations is 
obtained in section III. Since the spin orientation of particles in ordinary matter is random, the macroscopic spacetime 
average of the spin vanishes. In this case, the solutions of the cosmological equations are divided into two classes. Each 
of them is related with only one torsion function, the scalar or the pseudoscalar torsion function. They are obtained 
in section IV and V, separately, using different methods. For the scalar torsion function the equations take the form 
of a dynamical system, of which asymptotically stable critical points represent the exact de Sitter solutions. For the 
pseudoscalar torsion function an exact analytic solution of the cosmological equations is presented in section V. In 
terms of this solution the acceleration and the phase transformation from decelerating to accelerating expansion of the 
universe can be explained. All of these solutions indicate that in vacuum the spacetime possesses an intrinsic torsion 
which does not originate from the spin of matter. It is the torsion that causes the acceleration of the cosmological 
expansion in vacuum. The torsion of the spacetime can be produced by the energy and pressure besides the spin of 
matter. In section VI we obtain some exact analytic solutions of the cosmological equations in the case 7 = 0. These 
solutions can only describe the inflation (in the early epoch) or the decelerating expansion (in the later epoch) of the 
universe. This means that the term jT^ vp T is necessary to construct a model of cosmic acceleration. The section 
VII is devoted to conclusions. 



II. Gravitational field equations 



We start from the action 



-R + aR 2 + m^R^ + lT^ p T^P 



+ s m [sr^P^V] , 



(1) 



where I = U KG / c 3 is the Planck length, a, and (3 are two parameters with the dimension of I 2 , 7 is a parameter of 
dimensionless, ip denotes matter fields. In contrast with [27], here the connection is not symmetric, i.e. r a ^ 7 7^ r Q 7 ^, 
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and appears in the action of matter S m . In other words, we are dealing with a metric-affine theory rather than a 
Palatini one. By the same way used in [29, 30, 31], the variational principle yields the field equations for the tetrad 
ei p and the spin connection T IJ ^: 

& Rvp, ^ e ^R 

= E\-a {±e Iv R vp - e 1 p R) R-(3 (2e Ia ' R p ' a R pp + 2e J p RP° R 1 j ^ - e 1 \ l R pa R'" T ) 

+7 (4ft, (e IX T pX v ) - Ae K T e IX T^d v e K T + e\T x pa T x p ° - Ae Iv T x VT T Xp T ) , (2) 

e { i v ejfe K T d v e K T + e {I v e A T Y\ T + e {I »e,jfT x Xv 
= «„" - 4a (e^ejfT^rR + e^erf (T x Xv R - d v R) + e {I v ejfRe K T d v e K r ) 
-4/3e/ {e^d v R x v] + e^R^e K T d v e K T + ej^\ T R x ^ + e^R T ^Y% x ) 
-A ieiu ej T T^ T . (3) 

where E 1 ^ and sjj p are energy- momentum and spin tensors of the matter source, respectively. We use the 
Greek alphabet (/i, v, p 1 ... = 0, 1,2,3) to denote (holonomic) indices related to spacetime, and the Latin alphabet 
(/, J, K, ... = 0, 1, 2, 3) to denote algebraic (anholonomic) indices, which are raised and lowered with the Minkowski 
metric r\u = diag (—1, +1, +1, +1)- If a = (3 — 7 — 0, these equations become the field equations of Einstein-Cartan- 
Sciama-Kibble theory. Especially, (20) becomes the Einstein equation. To understand these equations, we will do 
a translation of (2, 3) into a certain effective Riemannian form-transcribing from quantities expressed in terms of 
the tetrad ej p and spin connection F IJ ^ into the ones expressed in terms of the metric g pll and torsion T x pll (or 
contortion K x pv ), as was done in [30]. It should be noted [16] that the set (e 1 p , T JK v ) corresponds to the first order 
formalism, while the set (<7 M „, T x pa ) to the second order formalism. The origin of this is that in the last case the 
non-torsional part of the affine connection is a function of the metric, while, within the gauge approach, the variables 
(e 1 M , T JK V ) are mutually independent completely. The descriptions in terms of the variables (e 7 M , T JK V ) and (g^, 
T x piJ ) are equivalent in our approach (the argument in detail see [16]). 

Subject to the metricity, the affine connection T X I1U is related to the tetrad ej p and the spin connection T IJ p by 

TV = e I x d p e I u + ej X e I ,r J Il , 

= (4) 
where { M A ;,}, K x pv are the Levi-Civita connection and the contortion, separately, with 

K x — —- (T x +T a +T x ) 

- 11 iiv — 2 V V v T ~ ± v[i ) i 

rp X p prpl _ pA _ pA 

T\ v = d p e\-dy P ^Y' j p e J Jv e\. (5) 
Accordingly the curvature R p apu can be represented as 

pp — p,pp J e 1 7 — ti v p - B v p + r*\ r A — v p \ r A 

— R p 4- 8 K p — hi K p +K P \ K x — K p \ K x 

+ {aM K x au - { x p v } K x ap + {A} K p Xp - { a x p } K p Xv , (6) 
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where R p ^ ailv = d^{ a p v } - d v {^ p ^,} + {a%} {<j X v} - {x p v}{a X ^} is the Riemann curvature of the Levi-Civita 
connection. In view of this, we can identify the actual degrees of freedom of the theory with the (independent) 
components of the metric g^ u and the tensor K x ^. 

III. Cosmological equations 

For the spatially fiat Friedmann-Robcrtson- Walker metric 

= diag (-1, a {tf , a (tf , a (i) 2 ) , (7) 
the non-vanishing components of the Levi-Civita connection are 

{o°o} = 0,{o ,} = {i° }=0,{i° j }=aaS ij , 

Wo} = 0,{/ } = {o i J -} = ^ ) {A}=0,i,j,fc,... = l,2,3. (8) 

The non-vanishing torsion components with holonomic indices are given by two functions, the scalar torsion h and 
the pseudoscalar torsion / [35] : 

Tiio = T220 = T330 = a 2 h, 

T123 = T 231 = T 312 = 2a 3 /, (9) 

and then the contortion components are 

"30 = 0, 
03 = h, 
22 = a 2 h, 
12 = -a/, 

-21 = 0/. (10) 
The non-vanishing components of the curvature R p ai iv and the Ricci curvature R^ v are 

#101 = R°202 = #303 = a 2 (h +H 2 + Hh+ h) , 

#123 = — #213 = #312 = 2a 3 / (H + h) , 

#203 = — #302 = #301 = -a ^Hf+ f^j , 

#212 = #313 - #323 - ({H + hf - / 2 ) , (11) 

Roo = -3 H -3 h -'AH 2 - 3Hh, 

Ru = R 2 2 = R33 = a 2 ^H + h+3H 2 + 5Hh + 2h 2 -f 2 y (12) 



#10 


= #20 


= #. 


#01 


= #02 


= #, 


IS® 
& 11 


= #22 


= K° 


#23 


= #31 


= K 3 


#32 


= #13 


= K 3 



R = 6H +6h +12H 2 + 18Hh + 6h 2 - 3/ 2 , 



(13) 
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where H —a (t) /a(t) is the Hubble parameter. Using these results and supposing the matter source is a fluid 
characterized by the energy density p, the pressure p and the spin s/j M we obtain four independent equations from 
(2) and (3): 

(H + h) 2 -f 2 - P - 

-(J3 + 3a) [4 (H + h) + »H (H + h) (h + h) 
-Ah (h + 2H) (h + H) 2 + 8(h + H) 2 f - 4/ 4 ] 

+2 7 (3h 2 + 4/ 2 ) = 0, (14) 

2 (h + h) + 3H 2 + AHh + h 2 -f-p 

+ {fi + 3a) [4 (H + h) + SH (H + h) (h + h) 
-Ah (h + 2H) (h + H) 2 + 8{h + H) 2 f + 4/ 4 ] 

-2 7 (2 h +8Hh + h 2 + 4/ 2 ) = 0, (15) 

(/3 + 6a) (n + h)+6{0 + 4a) {H + h) H + (5/3 + 18a) (H + h) h -A {(3 + 3a) / / 

+3 08 + 4a) hH 2 + (5/3 + 18a) h 2 H + 2 (/3 + 3a) h 3 - 2 {fi + 3a) fe/ 2 + ^h + ^soi 1 = 0, (16) 

f{2 (/3 + 6a) (h + h) + 6 (/3 + 4a) H 2 + 2 (5/3 + 18a) if/i 

+ (/? + 3a) (4/i 2 - 4/ 2 ) - 4 7 + 1} - l Sl2 3 = 0. (17) 

The system of the equations (14)-(17) has the similar structure as the system of gravitational equations for homoge- 
neous isotropic cosmological models in [31] except the coefficients. However, it is the differences in coefficients that 
make the system of the equations (14)-(17) easy to handle and possible to obtain some exact or analytic solutions in 
several cases as will be shown in the next sections. 
The equations (14) and (15) can be written as 

H + h=2 1 h -2H 2 + (87 -S)Hh- (2 7 + 1) h 2 + f + 1 (p + 3p) , (18) 

and 

(fi + 3a) [-4 (h + h) -8( y H + h)H(H + h) 
+Ah (h + 2H) (h + H) 2 - 8 (h + H) 2 f + 4/ 4 ] 

+ 7 (6^ 2 + 8/ 2 ) + H 2 + 2Hh + h 2 - f -\p = 0. (19) 

Since the spin orientation of particles in ordinary matter is random, the macroscopic spacetime average of the spin 
vanishes, we suppose s/,/ A = 0, henceforth. Then, the equations (16), (17) become 

(/3 + 6a) (h + h) + 6 (fi + 4a) (H + h) H + (5/3 + 18a) (H + h) h -A (/3 + 3a) / / 

+3 (fi + 4a) hH 2 + (5/3 + 18a) h 2 H + 2 (fi + 3a) h 3 - 2 (fi + 3a) hf 2 + ^h = 0, (20) 



and 



2 f{2 (J3 + 6a) (H + h)+6(P + 4a) H 2 



+2 (5/3 + 18a) Hh+(f3 + 3a) (Ah 2 - 4/ 2 ) - 4 7 + 1} = 0. 



(21) has the solutions 

/ = 0, 

and 



2 (/3 + 6a) /• -\ 3(/3 + 4a) (5/3 + 18a) 2 

f ~ WT3a) [ H + V + WT3a) H + WT3a) Hk + k 



(/9 +3a) 8(/3 + 3a)' 

We will solve the equations (18-20) in the cases (22) and (23), respectively in the next two sections. 

IV. Exact de Sitter solutions with scalar torsion function 

In the case / = 0, (18) and (19) can be written as 

1 



H= (27 - 1) h -2H 2 + (87 -3)Hh- (2 7 + 1) h 2 + ^ (p + 3p) . 



h + {4H _ h)h _L {H + hf , P + 3pY l2^ + 3a)(H + hr + 3(H + h) 2 + l 8h 2 ,- P 



and 



2 7 v ; 12 7 J 48 7 2 (J3 + 3a) 

which have the solutions 



■ 1 2 P + 3p / l2(/3 + 3a) (H + h) 4 + 3(H + h) 2 + l8h^-p 
h--(AH-h)h+-(H + h) "^-±y 48 7 2 (0 + 3a) ' 

H = -±(H + h) 2 -H 2 + 3Hh-h 2 + ^P 
27 127 

/ 12 (/3 + 3a) (ff + /i) 4 + 3 (H + /i) 2 + 18/i 2 7 - p 
±(27_1) V 48 7 2 (/? + 3a) ' 

Differentiating (24) gives 

H= (2 7 -l)h -AH H + (87 - 3) H h + (87 - 3) H h -2 (2 7 + l)hh+^(p +3 j>) . 
The equation (20) has the form 

(0 + 6a) (# + ft) + 6 (/3 + 4a) (H + h) H + (5/3 + 18a) (H + h) h 
+3 (/3 + 4a) hH 2 + (5/3 + 18a) h 2 H + 2 (/3 + 3a) /i 3 + = 0. 
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(28) and (29) have the solutions 



(48a 7 + 12/3-7 ~ 2/3) H + (4/3 7 - 3/3 - 6a) fe • 
17 ~ 2 7 (/3 + 6a) ff 

(2/37 - 2/3 - 12a 7 ) # + (14/3 7 + 60a 7 - 3/3 - 6a) /i 
2 7 (/3 + 6a) 

3 (2 7 - 1) (J3 + Aa)_ hH2 _ (2 7 -l)(5/3 + 18a) h2g _ 



2 7 (/3 + 6a) 2 7 (/3 + 6a) 

(2 7 -!)(/? + 3a) 3 _ 27-I J_ , ■ A 

^/3 + 6a) " 8^f/3 + 6ar + 12^V P ^J' W 



2/3H + (8 7 /3 + 48 7 a + 3,3 + 6a) h ■ 
k ~ 2 7 (/3 + 6a) H 

(2/3 + 8 7 /3 + 48 7 a) i? + (3/3 + 6a - 4 7 /3 - 24 7 a) h 



2 7 (/3 + 6a) 

3 (/3 + 4a) h 2 _ 5/3 + 18a 2 /3 + 3a 3 

" 2 7 (/3 + 6a) 2 7 (/3 + 6a) 7 (/3 + 6a) 

■wW"^( p_3p )' (31) 



Letting 



H= X, h= Y, (32) 

we have the dynamical system 

• _ (48a 7 + 12/37 -2/3)H+ (4fi 7 - 3/3 - 6a) h 
X ~ 2 7 (/3 + 6a) 

(2^ 7 - 2/3 - 12a 7 ) g + (14/3 7 + 60a 7 -3/3-6a)h 

2 7 (/3 + 6a) y(ii '^ 
3 (2 7 - 1) (/3 + 4a)_^ ff2 _ (2 7 -l)(5/3 + 18a) ^ 2g _ 



2 7 (/3 + 6a) 2 7 (/3 + 6a) 

] -h*-*-\ > + 



(2 7 -l)(/3 + 3a) _ 27-1 J_ / • A 

7 (/3 + 6a) 8 7 (0 + 6a) 12 7 V J ' 



ft = y (i?, h) = -^-(H + h) 2 - (AH — h)h— y^- (p - 3p) 



12(1 + 3a) (g + fef + 3 (H + fe) 2 + 18 7 fe 2 - p 
V 48 7 2 (/3 + 3a) ' W 



The critical point equations consist of 



-3 (/3 + 4a) Mf 2 - (5/3 + 18a) ft 2 i? - 

= o, 



X = 0, 
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(H + ft) 2 - 2 7 (AH - ft) ft - 1 (p - 3p) 

/ i2(/3 + 3a) (ff + ft) 4 + 3(g + ft) 2 + 18 7 ft 2 -p _ n 
± V 12(0 + 3a) (d4) 

In order to discuss the stability of the critical points we need to calculate the matrix elements of the Jacobian: 

d X (48a 7 + 12/3 7 -2/3)H+ (407 - 3/3 - 6a) ft 

ax ~~ 2 7 (/3 + 6a) ' 

dX _ (48a 7 + 12/3 7 - 2/3) (2/3 7 - 2/? - 12a 7 ) 

2 7 (/3 + 6a) 2 7 (/3 + 6a) 1 ' j 

(2/3 7 - 2/3 - 12a 7 ) if + (14/37 + 60a7 - 3/3 - 6a) ft dY (H, ft) 

2 7 (^ + 6a) a# 
3 (2 7 - 1) (/3 + 4a) (27 - 1) (5/3 + 18a) 2 

7 (/3 + 6a) 2 7 (/3 + 6a) ' 



91 = (4/3 7 - 3/3 - 6a) (14/3 7 + 60a 7 - 3/3 - 6a) 
9ft " 2 7 (0 + 6a) 2 7 (/3 + 6a) [ ' j 

_ (20 7 - 2/3 - 12a 7 ) if + (14/37 + 60a7 - 3/3 - 6a) ft dY (H, ft) 
2 7 (0 + 6a) 9ft 
3(2 7 -l)(/3 + 4a) g2 _ (2 7 -l) (5/3 + 18a) 

2 7 (/3 + 6a) 7 (/3 + 6a) 

3(27- 1) 08 + 3a) 2 2 7 - 1 



7 (/3 + 6a) 87 (/3 + 6a) 



^ = 1^=0^ = 0^=0 

ax ' an ' dh ' ax ' 



aft _ ^(^ ) = i /i_ 



a# aF 7 V7 

8 (/3 + 3a) (H + hf + H + ft 



4 7 ^4 (/3 + 3a) 2 (ff + ft) 4 + (/3 + 3a) (if + ft) 2 + 6 (/3 + 3a) 7 ft 2 - ± (/3 + 3a) p' 







(H* + 


(H 







aft aft 

8 Qg + 3a) (if + ft) 3 + H + ft + 6 7 ft 

± (35) 

4 7 y 4 (/3 + 3a) 2 (if + ft) 4 + (13 + 3a) (if + ft) 2 + 6 (/3 + 3a) 7 ft 2 - ± (/3 + 3a) p 

In order to stress the role of the torsion as the source of the accelerating expansion of the universe we concentrate 
on the vacuum solutions for some special choices of the parameters a, /3, and 7. The critical point equations (34) can 
be simplified and solved exactly when /3 — 4a or /3 = 3a. 
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A. When /3 = -4a 

In this case the gravitational Lagrangian is a special case of quadratic curvature gravities [36] when the torsion 
vanishes. 

In vacuum the dynamical system (33) becomes 

X = \4h--H- J-ft | X+[5H-h--H- J-ft ) Y (H, h) 
V 7 2 7 ) \ 7 2 7 ) 

, (27"1) ; 



167a 



-h(8hHa + 8ah 2 - l) 



H= X, 



h = Y(H,h) = ^-(H + hf - (AH - ft) ft 
27 

± 1 Ua(H + hf -(H + hf -'^jh 2 
47 V a 



and the critical point equations (34) become 



ft (8hHa + 8ah 2 - l) = 0, 



X = 0, 



(37) 
(38) 



The equation (37) 



leads to 



(H + hf -2 7 (4H-h)h 



1 4a(H + hf -(H + hf -6jh 2 _ Q 
2V a 



ft (8ft.Ha + 8aft 2 - l) = 0, 



ft = 0, 



(39) 



or 



Then we have two cases. 

In the first case the solution 



8hHa + 8ah 2 -1 = 0. 



h = 0,H = 0,H= X = 0, 



(40) 



corresponds to a static Minkowski spacetimc. 
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In the second case, h and H satisfy the equations 



8hHa + 8ah 2 -1 = 0, 

and 

(H + hf - 2 7 {AH - h) h 



1 Ua(H + hf - {H + hf-^K 2 n 
± 2 V a = °' 



which can be written as 

1 



H = - h+ Sah> 



and 



25600of 3 /iV + 128 7 a 2 /i 4 (3 - 40 7 ) + 16jah 2 (I67 + 5) - 87 + 1 = 0. 

In order to obtain a concrete results we give some specific value of 7. 
When 



(41) 



7 = 4, 



(42) 



the equations (41) become 



and 



H = -h + 



8aft' 



409600a 3 /i 6 - 80384aV + 4416a/i 2 - 31 = 0, 



which have a real solution 



TT 1.40162 , 0.0841322 
H = = — , h = 



The Jacobian matrix of the dynamical system (36) given by (35) is 



M 



I 0. 39583 2.093 33. 382 \ 

^/a a a 

1 



33904 5.4102 i 



which has the eigenvalues: — 0. 85371/^/a, — 4. 951/y/a, — 1. 3522 x 10 3 /y/a. If a > 0, all the real parts of the 
eigenvalues are negative. This means that the critical point 



n „ 1.40162 0.0841322 
X c — 0, H c — -= — , n c — 



y/a 



(43) 
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is asymptotically stable and then 

H= X = 0, (44) 

gives an asymptotically stable de Sitter solution. 
By the same way we can compute for 



7 = 2, (45) 

the dynamical system (36) has a real critical point 

„ TT 0.824837 , 0.130802 
X c = 0,H c = j^,h c = (46) 

the Jacobian has the eigenvalues: -1. 665/Va + 0. 26071i/v / a, -1. 665/Va 7 - 0. 26071i/V<*, -3. 6641 x 10~ 3 /\/a- 
For 

7=1, (47) 

(36) has a real critical point 

„ TT 0.488003 , 0.185576 , . 

X c = 0,H c = 7 ^,h c = 48 

the Jacobian has the eigenvalues: -0. 8932/Va + 0. 97558i/V^ -0. 8932/V« - 0. 97558i/Va, -0.04017/^. 
For 

1=\ (49) 



(36) has a real critical point 



the Jacobian has the eigenvalues: —2. 6618/ \/a, —0. 22666/-y/a, —3. 5927 x 10 5 /\/a. All the three critical points are 
asymptotically stable. 
For 

7= \, (51) 



(36) has a real critical point 



0.488003 , 0.185576 , N 

X c = 7 H c = y^,h= 52 



the Jacobian has the eigenvalues: —5. 1634/y / a, —0. 68424/- s /a, 3. 8689 x 10~ 2 / \fa. One of the eigenvalues is positive, 
which means that the critical point (52) is unstable. These examples illustrate that the stability of the critical points 
depends on 7, i.e. on the term r yT^ l/p T^ l/p in the action (1). 
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B. When /3 = 3a 

This corresponds to conformal (Weyl) gravity (resent see [37]) and Critical Gravity [38] when the torsion vanishes. 
In this case the system of the equations (14)-(17) has the form 

3H 2 + 6Hh + 3 (67 + 1) h 2 + 3 (87 - 1) f - p = 0, (53) 
2 H +2 (1 - 2 7 ) h +3H 2 + 4(1- 4 7 ) Hh+{1- 2 7 ) h 2 - (1 + 87) f - P = 0, (54) 
3a (h + hj +6a(H + h) H +3a (H + h) h +3aHh (H + h) + jh+ ^soi 1 = 0, (55) 

/{6a (H + h)+ QaH 2 + 6aHh - 4 7 + 1} - ^s 12 3 = 0. (56) 
When su^ = 0, the equation (56) leads to 

/ = 0, 

6a (Jl + h)+ QaH 2 + QaHh - 4 7 + i = 0. (57) 



We deal with only the first case / = in this section, then the equations (53)-(55) can be written as 

-H± J-6~{H 2 + i (I + 67) p 

h = TT~a ' 

I + 67 



(58) 



h= H +^ tH 2 + ?—^Hh - -h 2 - -— 1 -p, (59) 

2 7 -l 2(2 7 -l) 27-I 2 2(27- 1)^' v ' 

and 

H + h+2(H + h)H + (H + h)h +hH (H + h) + -^-h = 0, (60) 

12a 

Differentiating (59) gives 

1 3 ■ 2-87 ■ 2-87 ■ ■ 1 

h =—i H+ — HH+ —i hH+ ^ Hh - hh -WTT) P - (61) 



Substituting (58), (59) and (61) into (60) yields 
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where 



Letting 



H 



2 (12 7 2 - 87 - 3) 
' (2 7 -1) (I + 67) 
(636 7 2 - 48 7 - 7 



H H + 



(67 +1) 



P H 



rH 6 



(67-I) (26 7 - 5) „ 2 



H Z P- 



IO7-3 
47(1 + 67)" 



2 (2 7 - 1) (1 + 67)" (2 7 - i)(i+6 7 r 

(2 7 -l)g (2 7 -l)P 67-I 1 ■ 

24 7 a(67+l) 24 7 a(67+l) 4 7 (2 7 -l) 1 4 7 ' 



T HP A 



p = ± x l~& 1 m + -{i + Q 1 ) P . 



(62) 



(63) 



H= X, 



(64) 



we have the dynamical system 



H = X, 

■ = 2(l2 7 2 -8 7 -3) 2 

(2 7 -1) (67 + I) (67 + I) 



FX 



(636 7 2 - 48 7 - 7) 



ff3 (67 - 1) (26 7 - 5) 



PH Z ~ 



IO7-3 
47(67 + 1)' 



2(2 7 -l)(6 7 +ir (2 7 -l)(6 7 +l) 2 

(2 7 -l)g (2 7 -l)P 67-I l_ p 

24 7 a (67 + 1) 24 7 a (67 + 1) 4 7 (2 7 - 1) P 4 7 ' 



r P 2 ff 



with the matrix elements of its Jacobian: 



(65) 



dH 
dH 

dX 
dX 

dX 
dH 



' dX ' 
2 (12 7 2 - 87 - 3) 
" (2 7 -1) (I + 67) 
2 (12 7 2 - 87 - 3) 



# + 



X 



(67 + I) ' 
12 7 i? 



X 



(2 7 -l) (I + 67) (6 7 +l)P 
67(67-1) (26 7 - 5) ff 3 + 3 (676 7 2 - 8O7 - 1) g2 
(2 7 -l)(l + 6 7 ) 2 P 2(2 7 -l)(l + 6 7 ) 2 



2 (67 - 1) (26 7 - 5) 
" (2 7 -l)(l + 6 7 ) 2 



HP 



IO7- 3 
4 7 (1 + 67)" 



rP A 



+ 



2 7 



+ 



2 7 -l H 



67-I 



247a (67+I) 4a(67+l)P 47(27-1) 



(66) 



The fixed point equations are 



X = 0, 

(636 7 2 -48 7 -7) g3 _ (6 7 -l)(26 7 - 5) pg2 _ 10 7 - 3 p2R 
2 (2 7 - 1) (67+ l) 2 (2 7 -l)(6 7 +l) 2 4 7 (6 7 + l) 2 

I (27~l)g (27-1)^ 67-I 1 ■ 

24 7 a(67+l) 24 7 a(67+l) 4 7 (2 7 -l) ^ 4 7 
= 0. (67) 



Using (63) we obtain the equation 



348 7 2 -48 7 +l |g3 | (2 7 ~l)g ( (6 7 -l)(26 7 - 5) g2 , (2 7 -l) 



(2 7 -l)(l + 6 7 ) 2 24 7 a(6 7 +l) ^ (2 7 - 1) (6 7 + l) 2 24 7 a(6 7 +l)^ 

10 7 -3 „ fry-1 „ 1 ■ 



12 7 (6 7 + 1) 4 7 (2 7 -l) 4 7 

= 0. 



In vacuum 



p = p = 0, P = vH, v = ±y/^6j, 

The fixed point equation (68) becomes 

348 7 2 - 48 7 + 1 g3 (2 7 -l)ff _ / (6 7 -l)(26 7 - 5) g2 2 7 - 1 \ ^ 

(2 7 - l)(l + 6 7 ) 2 24 7 a(6 7 +l) ^ (2 7 - 1) (6 7 + l) 2 24 7 a(6 7 + l)J 

which leads to 

H = 0, 



or 

H 2 = (^-l)(2 7 -l) 2 (l + 6 7 ) 

24 7 a ((156 7 2 - 56 7 + 5) v - 348 7 2 + 48 7 - 1) ' 

In the first cas, we have 

H= X = 0, H = 0, h = 0, 

which correspond to a static Minkowski solution. 
In the second case 

H 2 = __ (^-l)(2 7 -l) 2 (l + 6 7 ) 



24 7 a ((156 7 2 - 56 7 + 5) v - 348 7 2 + 48 7 - 1) ' 
according to (66) the Jacobian matrix has the form 



M 



1 

b c 



with 



b = 6 7 (6 7 ~l)(26 7 - 5) ff 2 | 3 (676 7 2 - 80 7 - 1) g2 
(2 7 - l)(l + 6 7 ) 2 v 2(2 7 -l)(l + 6 7 ) 2 
2(6 7 -l)(26 7 - 5) ^ g2 _ 10 7 -3 i 2r2 
(2 7 -l)(l + 6 7 ) 2 4 7 (l + 6 7 ) 2 

2 7 - 1 2 7 - 1 



24 7 a(6 7 + l) 4a(6 7 + l)z/ 
2 (12 7 2 - 8 7 - 3 - 2z/ 7 + v) 

(2 7 -l)(l + 6 7 ) H ' 

and has the eigenvalues: \c + |y / (c 2 + 4b), \c— \ \J (& + 46). 
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If 



(69), (74) and (76) give, separately, 



7 24' 



V 2' 



H 2_ 169 
948a' 



408901 _ 13 / 169 
~~ ~ 30336a ' C ~~ ~TV 948a' 

Then the Jacobian matrix has the eigenvalues: 



1 1 n~o 7TT 13 /169 i 19807661 

2 C+ 2^ (C +46) = -lV^VM8 + V^V"^7iW' 

Ic-l^TW) 13 m 1 1 19807661 



2 2 VV ' 14 V / SV 948 2 V / ^V 371616 

If a. > 0, all the real parts of the two eigenvalues are negative. This means that the critical point 



169 

X c = 0,H c = \l— — , (77) 
V 948a v ' 

is asymptotically stable and then 

H= X = 0, (78) 

gives an asymptotically stable de Sitter solution. (58) gives 

13 ^237 
711 

oi- 



lS V237 

h = =-. 

237 V" 

We have seen that by appropriate choices of 7, we can obtain asymptotically stable de Sitter solutions in both 
cases, /3 = 4a and (3 = 3a, though the cosmological equations have different structure. This means that the structure 
of the cosmological equations depends on /3, while the stability of the solutions depends on 7. 

In all of the solutions obtained the torsion function h does not vanish in vacuum, which means that the torsion is 
an intrinsic geometric nature of the spacetime. It is the torsion that causes the accelerating expansion of the universe 
in vacuum. 
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V. Analytic solutions with pseudoscalar torsion function 

Differentiating (23) gives 

■ /3 + 6a /•• -\ 3(/3 + 4a) rr • 5/3 + 18a • 5/3 + 18a • , • 

ff=T73 ST # + ft + 7773 # +777 rr # ^ + 77-3 77^ ft +ft ft . 79 

4(/3 + 3a) V / 2(/3 + 3a) 4 (/3 + 3a) 4 (/3 + 3a) 

Substituting (23) and (79) into (20) gives 

ft = 0. (80) 

Then the equations (18), (19) and (23) become 

H=-2H 2 + f+ l -{p + Zp), (81) 

(/3 + 3a) [-4 H -8HH 2 - 8H 2 f + 4/ 4 ] + (87 - 1) f + H 2 - \p = 0, (82) 
f2 _ (ff + 6a) r > , 3(/3 + 4a) 7 1 

j 2(/3 +3a) 2(/3 +3a) (/3 + 3a) 8 (/3 + 3a) ' k j 

They have the solutions 



H _ _ gl^M _ _L p + 14%- „ + 3p) _ 0i±M|±M „ + ^ (8 , 



2 I-87 1 167(/3 + 3a)-3(/3 + 4a) (/3 + 3a) (/3 + 4a) 2 
J = W + ^ 48^3 (P + 3P) + 24^ (P + 3P) ' (86) 

The equations (84) and (85) ply the roles of the Friedmann equation and the Raychaudhuri equation in General 
Relativity The equation (86) indicates that even in vacuum the spacetime possesses the torsion / = yj X 3 ~^J , which 
has been found in [32] . Hence the conception of the vacuum as physical notion is changed essentially. Instead of it 
as passive receptacle of physical objects and processes, the vacuum assumes a dynamical properties as a gravitating 
object. The combination of (84) and (85) yields the acceleration equation 



a 87-I /3 +3a 

a = — W + ^T iP + 3p) - (8?) 

Letting 



/3 = na, 
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wc have 



a 1 - 87 n + 3 . 

+ (P + 3p) . 



a Ana 3n 
Some important consequences can be obtained from (88): 

i) The term 1 4 ^ y plies the role of the cosmological constant, which agrees with the result in [31]. If 1 4 ~^ 7 > 0, 
p = p = 0, then a> 0, the acceleration of cosmological expansion acquires the vacuum origin. 

ii) If 

n > 0, or n < -3 (89) 
p + 3p accelerates the expansion of the universe. If 

- 3 < n < 0, (90) 

p + 3p decelerates the expansion of the universe. Especially, when n = —2,7 = 1/8, (88) becomes the acceleration 
equation in general relativity. In other words, the latter is only a special case of the former. 

iii) If 

n > 0,7 > 1 or n < -3, 7 < i, (91) 

the universe can undergo a phase transformation from an accelerating to a decelerating expansion. 

iv) If 

-3<n<0, 7 >i (92) 

the universe can undergo a phase transformation from a decelerating to an accelerating expansion. 

We find this picture very appealing and physical since it seems to indicate that in metric-affine gravity as matter 
tells spacetime how to curve, matter will also tell spacetime how to twirl. 

VI. Analytic solutions in the case 7 = 

In the last two sections we have seen that the coefficient 7 plies a central role in determining the behavior of the 
scale factor and the evolution of the universe. In order to investigate this point thoroughly, we discuss a extreme case, 
7 = 0. In this case the equations (14)-(17) take the form 

(H + h) 2 -f 2 :; 
+ (0 + 3a) [-4 (h + h) - SH (H + K)(h + h) 

+4h (h + 2H) (h + H) 2 - 8 (h + H) 2 f + 4f 4 } = 0, (93) 



2 (h + h) + 3H 2 + AHh + h 2 - f-p 

2 

-{/3 + 3a) [-4 (H + h) - 8 (H 2 + Hh) (h + h 



+Ah (h + 2H) (h + H) 2 -8{h + H) 2 f + 4f 4 } = 0, (94) 
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08 + 6a) (h + ft) + 6 (/3 + 4a) (H + ft) H + (5/3 + 18a) (H + ft) ft -4 (/3 + 3a) / / 

+3 (P + 4a) hH 2 + (5/3 + 18a) h 2 H + 2 (/3 + 3a) ft 3 - 2 (/3 + 3a) ft/ 2 + -ft + ^.soi 1 = 0, (95) 

f{2 {p + 6a) (h + ft) + 6 {p + 4a) H 2 

+2 (5/3 + 18a) Hh+{p + 3a) (4ft 2 - 4/ 2 ) + \} - \si2 S = 0. (96) 



(93) and (94) can be written as 



H + h= -2H 2 - 3Hh - ft 2 + f + \{p + 3p) , (97) 



3 (ft + H) 2 - 3f - p 

+4 (p + 3a) ((if + ft) 2 - / 2 ) (p + 3p) - i (/3 + 3a) (p + 3pf 
= 0, (98) 

Differentiating (97) gives 

H + h= -4H H -3ft H -3H ft -2ft ft +2/ / + 1 (p +3 p) . (99) 
Substituting (97) and (99) into (95) and (96) yields 

-P (ft + 4£T) (ft + H) 2 + p (ft + 2ff) / 2 - 2/3/ / +h + 1 (2/3/i + 6aft + 3/3ft) (p + 3p) 

+ i(/3 + 6a) (p +3p) + \ s oi 
= 0. (100) 



and 



/{2/3 (ft + if) 2 - 2pf + \(P + 6a) (p + 3p) + 1} - i Sl2 3 = 0. (101) 



If 

= 0, 

(101) leads to 

/ = 0, 

or 



2/3 (ft + H f - 2/3/ 2 + 1 (p + 6a) (p + 3p) + ^ = 0. 



(102) 



A. When / = 



The equations (98) and (100) become 

(H + h) 



■ 2 + 3a) (p + Spf + p 

4(/3 + 3a)(3p + p) + 3 ' 

and 

1, 1 



which have the solutions 
H 



■gh+- (2/3H + 3f3h + 6ah) (p + 3p) + (/3 + 6a) (p +3 p) 
P(h + 4H) (H + h) 2 = 0, 



f + f g (p + 5p) + 6a (p + 3p) + \(P + 3a) (5,3 + 12a) (p + 3p) 2 



| + |/3 (p + p) + 6a (p + 3p) + i (/3 + 3a) (5/3 + 12a) (p + 3p) 2 
+ 



4(/3 + 3a) (g + 6a) (p + 3p) + 3 Q8 + 6a) 

3 i 9 a ( „ i „\ i f „ i Q,^ i I I a _i_ 9„A /c a i i 9_,\ ^ _i_ l^ 2 



, 2 /S (p + P) + 6a (p + 3p) + i (/3 + 3a) (5/3 + 12a) (p + 3p) 
3/3 (p-p) 



(p +3 p) , 



| + p (p + p) + 6a (p + 3p) + § 08 + 3a) (5/3 + 12a) (p + 3p) 2 

4 (/3 + 3a) (/3 + 6a) (p + 3p) + 3 (g + 6a) 

| + |/3 (p + p) + 6a (p + 3p) + ± (/3 + 3a) (5/3 + 12a) (p + 3p) 



J 

(p +3 p) 



with 



§Q3 + 3a ) (P + 3 P) 2 + P 
4 (/3 + 3a) (3p + p) + 3 ' 



J = ±\ 

The equation (97) now becomes 

H + h=-(2H + h) (H + h) + l(p + 3p) . 

o 

Letting 

/3 = na,p = wp, 

(105), (106) and (107) can be written as 



_ (f + f^ap + \Ba 2 p 2 ) J + (ADap + G) a (1 + 3w) P 
H= | + |Cap+iBa2p2 ' 



3n (1 - w) apj - (4L>ap + G) a (1 + 3w) P 
f + |Gap+iBaV ' 



J = ± V^T3 /5 ' 
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whore; 

A = n + 4 + (5n + 12) w, B = (n + 3) (5n + 12) (1 + 3w) 2 , 
C = 3n + 4 + (3n + 12) to, D = (n + 3) (n + 6) (1 + 3to) , 

E = (n + 3) (1 + 3w) 2 , F = (n + 3) (3w + 1) , G = 3 (n + 6) . (113) 
The equation (108) gives the acceleration equation 



- = - h-(H + hf -hH +l(l + 3w)p, (114) 
a o 



Let us consider two special cases, 
i) For the early universe, 

ap » 1, 

we compute using (110)- (114) and obtain approximately 



a 12 

This represents an inflation universe, 
ii) For the later epoch 



" i(l + 3«0p. (H5) 



wc; have 



ap <C 1, 



^ = -ip(l-3t«). (116) 



This represents a uniformly expanding universe if w = 1/3 (radiation epoch), or an decelerating universe if w < 1/3 
(matter epoch). 



B. when / =£ 



The function / satisfies the equation (102), which yields 



f 2 = (H + h) 2 + ^(p + 3 P ) + ± (117) 



and 



Using (97) we have 



H + h=-H(H + h) + P±^(p + 3p) + ± (118) 
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ff=-H(H + hf + ^±^(H + h)(p + 3p) + ^(H + h) + ^^(p+3p) 1 (119) 
substituting in to (98) and (100) we have 

-^T (P + 3P) ~ WWT 3a) iP + 3P) ~ WT3aj 4/3 {? + 3a) = °' (120) 

= 0. (121) 

So the field equations have no definite solution. 

The results obtain above indicates that if 7 = in both cases, / = 0, and / ^ 0, there exists no solution describing 
an accelerating universe. In other words, the term ')T IJ 'vpT ll vp is necessary to the existence of the solutions describing 
an accelerating universe. 



VII. Conclusions 



Quadratic theories of gravity described by the Lagrangian R + aR 2 + l3Rp V R pl ' have been studied in many works 
in supergravity, quantum gravity, string theory and M-theory. However, the cosmology in these theories has not been 
explored extensively, especially, when the torsion of the spacetimc is considered. In this paper we show that by only 
allowing the connection to be asymmetrical and adding a term ^T p vp T p vp to the Lagrangian R + aR 2 + f3R tlu R pl ' 
some meaningful cosmological solutions can be obtained. These solutions provide several possible explanations to the 
acceleration of the cosmological expansion without a cosmological constant or dark energy. One can find that although 
the field equation (2) returns to Einstein's equation when a — (3 = 7 = 0, the cosmological equations (18-21) are 
essentially different from the Friedmann equation and the Raychaudhuri equation and then give different description 
to the evolution of the universe. The acceleration equation of the universe in general relativity is only a special case 
of the equation (87). These equations involving higher-derivatives can be solved by appropriate choice of /? and 7. 
Not only numbers of asymptotically stable de Sitter solutions expressed by critical points of a dynamical system but 
also exact analytic solutions are obtained. These solutions indicate that the terms f3R lxv W lv and jT p up T ' ^ p ply the 
different roles: the former determines the structure of the equations while the latter determines the behavior and 
the stability of the solutions. To construct a model of cosmic acceleration the Lagrangian R + aR 2 + f3R lllJ R p,u + 
"fT^vpT^P is sufficient and necessary. 

Owing to the solutions obtained some conceptions have to be changed essentially. According to these solutions, 
even in vacuum the spacetime can possesses torsion and curvature. Therefore, instead of vacuum as passive receptacle 
of physical objects and processes, the vacuum assumes a dynamical property as a gravitating object. It is the torsion 
of the spacetime that causes the acceleration of the cosmological expansion in vacuum. Both the torsion and the 
accelerating expansion possess geometrical nature and do not invoke any matter origin. Furthermore, the energy and 
pressure of the ordinary matter can produce the torsion of the spacetime and cause either the deceleration or the 
acceleration of the cosmological expansion depending on choices of /3 and 7. 
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